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Triad interactions in the dissipation range 


By S. Kida 1 , R. H. Kraichnan 2 , R. S. Rogallo 3 , F. Waleffe 4 , and Y. Zhou 4 


Nonlocality of the triad interactions in the dissipation range of developed turbu- 
lence is investigated by numerical simulation and the quasi-normal theories. It is 
found that the energy transfer is dominated by nonlocal triad interactions over the 
wavenumber range 0.1 < k/k j < 4, where k$ is the Kolmogorov wavenumber. The 
nonlocality of the interaction has a close relation with the power of an algebraic 
prefactor of the exponential decay of the energy spectrum in the far-dissipation 
range. 


1. Introduction 

The triad interaction is the fundamental coupling among the various Fourier 
components of a turbulent velocity field and transfers energy predominantly from 
lower to higher wavenumbers. Properties of the triad interactions were studied by 
Kraichnan (1971) and have recently been analyzed numerically using data from 
numerical simulations (Domaradzki 1988; Yeung Brasseur 1991; Domaradzki & 
Rogallo 1990; Ohkitani & Kida 1992) and by an analysis of interactions among 
helical waves (Waleffe 1992). 

In the inertial range where the energy spectrum obeys the Kolmogorov —5/3 
power form, 

E(k) = C K e 2 ' 3 k~ 5/3 , (1.1) 

where e is the energy dissipation rate and Ck is the Kolmogorov constant (Kol- 
mogorov 1941), the interaction is local in the sense that triad interactions of scale 
disparity (see (3.4) for definition) greater than 10 (100) contribute only 15% (1%) 
of the energy flux (Kraichnan 1971; Ohkitani &. Kida 1992). Although the locality 
of the triad interaction is very weak, it is compatible with Kolmogorov’s (1941) 
concept of local energy cascade. 

The nature of triad interactions in the dissipation range is expected to be different 
from that in the inertial range because the energy spectrum decreases very rapidly 
(exponentially) in the dissipation range. In this paper, we investigate the interac- 
tions by analyzing numerical turbulence and by its comparison with the prediction 
of the quasi-normal closure theories. 
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2. Fundamentals 

We consider the energy dynamics of the motion of an incompressible viscous fluid 
in an infinite domain. The equation of motion is written in the Fourier representa- 
tion as 

«k(p)w?(q)-v* 2 u/(k) + £(k) (2.1) 

p+q+k=0 

with the continuity equation 

kj Sj( k) = 0, (2.2) 

where u,(k) is the x } (j = 1,2,3) component of the Fourier coefficient of velocity 
at wavevector k, 


JVmOO = h (s„ -*£■) + * (fik - tjzr) (2.3) 

is a third order tensor, v is the kinematic viscosity, and fj( k) is the external forcing. 
Here the time argument t is omitted for brevity, the asterisk denotes the complex 
conjugate, and repeated subscripts are summed over 1 to 3. 

The energy spectral density at wavevector k, 

E(k) = I|u(k)| 2 (2.4) 

evolves according to 

^E( k) = T(k) - 2uk 2 E(k) + Ro {uj(k)^(k)} , (2.5) 

which is derived by multiplying (2.1) by u*j and taking the real part. 

The first term on the r.h.s. of (2.5) represents the rate of energy transfer to the 
Fourier mode k through the nonlinear interactions with all the other modes, the 
second the energy dissipation by the viscosity, and the third the energy input by 
the external force. 

The energy transfer function T(k) is written as 

7\k) = £T(k|p,q), (2.6) 

p.q 

where 

T(k|p,q) = -ilm {P>*((k)u>(k)u*(p)£,(q)} £ k + P +q (2.7) 

is the triad energy transfer due to the interaction among three wavevectors k, p, 
and q that constitute a triangle (k + p + q = 0). Through a triad interaction, 
energy is exchanged among the three modes involved, with the total energy being 
conserved. That is, the following detailed balance of energy holds; 


r(k|p, q) + T(p|q, k) -I- T(q|k, p) = 0. 


( 2 . 8 ) 
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Consequently, the sum of the energy transfer function over all the wavevectors 
vanishes, 

£T(k)= ]T T(k|p,q) = 0. (2.9) 

k k,p,q 


Since the flow field is statistically isotropic, it is convenient to average each term in 
(2.5) over spherical cells in the wavevector space. We introduce the band-averaged 
energy spectrum 

E(k)Ak = Y E ( k ')’ ( 2 - 10 ) 

*-$Ai<|k'|<k+jAfc 


the band-averaged energy transfer function 


T(k)Ak = Y, T ( k ')’ (2* 11 ) 

k-$Ak<\k'\<k+\Ak 


and the band- averaged forcing spectrum 

F(k)Ak = Y Re{«*(k')/>(k')}. (2.12) 

*-}A*<|k'|<Jfc+$A* 


The energy spectrum equation (2.5) is then written for the band-averaged quantities 
as 

■^■E(k) = T(k) — 2vk 2 E(k) + F(k). (2.13) 

at 

The triad energy transfer function T(k|p, q) is also averaged over a spherical cell 
as 

f(k\p,q){Akf= Y r(kV,q')- (2-14) 

Jt — -j Ai<|k'|<fc+5 Afc 
p-\ Ak<\p'\<p+\Ak 
Afc<|q'|<*+$Afc 

The energy transfer function is then written as 

T(Jb) = (Afc) 2 ^f(*|p,$). (2.15) 

PiQ 


The detailed balance of energy (2.8) and the conservation of total energy by all the 
triad interactions (2.9) are written respectively as 

T(k\p,q) + T(p\q , k ) + T(q\k,p) = 0 (2.16) 

£f(fc) = °- 


and 


(2.17) 
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FIGURE 1 . Energy Spectrum in the dissipation range. The decay with wavenumber 
is essentially exponential. 

3. Numerical Simulation 

The equations of motion (2.1) and (2.2) are solved numerically in a periodic cube 
with an instability forcing, 

/>(k) = Cujik) for |k| < fc 0 . (3.1) 

The constant C is chosen to force the flow field to equilibrium with a specified 
range of forced scales ko = 3 and a specified range of resolved scales, in Kolmogorov 
units, fc m *xAd = 4. The Fourier- spectral method (Rogallo 1981) is employed for 
the spatial resolution, and time is advanced with a second-order Runge-Kutta pro- 
cedure. The alias errors arising in the nonlinear terms is removed by a combination 
of coordinate shifting and spectral truncation. The computational mesh (in physi- 
cal space) is 256 3 . The initial flow field is taken after over five large-scale turnover 
times from forced turbulence created at a lower resolution (128 3 ) at about the same 
Reynolds number (R\ % 65). The 256 3 field was then advanced until an equilibrium 
between transfer, and dissipation was achieved at the higher wavenumbers. We will 
study that transfer here in some detail. The Kolmogorov dissipation wavenumber, 

£d = (v 3 /e) lf4 (3.2) 

is about 30, so that the maximum resolved wavenumber Jt m „ = 121 retained in the 
simulation is about four times the Kolmogorov wavenumber. 

S.l Energy spectrum 

The energy spectrum at the final time of the simulation is shown in a semi- 
logarithmic plot in figure 1. The nearly straight line indicates that the energy 
spectrum decays essentially exponentially with wavenumber. In order to examine 
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FIGURE 2. Logarithmic derivative of the energy spectrum function. 

For E(k ) = A(k/k d ) 0 exp[-l3k/k d } the slope of the curve gives -/? and the in- 
tersection on the vertical axis gives a. , simulation data; , least-square 

fit over .5 < k/k d < 3 gives a = —1.6, /? = 5.2. 

the shape of the spectrum more precisely, we compare it with an exponential form 
with an algebraic prefactor 

E(k ) = A(fc/fca)“exp[— (dk/ki]. (3.3) 

In figure 2, we plot the logarithmic derivative of the energy spectrum. If the spec- 
trum has the form (3.3), the slope of the line and its intersection with the verti- 
cal axis give the values of 0 and a, respectively. Since this is an instantaneous 
spectrum, the fluctuations are quite large. Nevertheless, the least squares fit over 
0.5 < k/k d < 3 gives a = -1.6 and /? = 5.2, which are consistent with those 
observed before in numerical simulations by other researchers (Kida & Murakami 
1987; Kida et al. 1990; Kerr 1990; Sanada 1992). More recently however, Chen 
(1992 private communication) has simulated the dissipation range at R\ ~ 15 to 
much higher k/k d . He finds a = 2.16 and /? = 7.35 by a least-square fit over the 
range 5.2 < k/k d < 10.4. The data from the present simulations do not coincide 
with Chen’s for k/k d < 4, suggesting that the results are sensitive to Reynolds 
number. There may also be some effect due to the method of forcing. Incidentally, 
the exponential shape of the energy spectrum in the far-dissipation range has also 
been observed in experiments (Sreenivasan 1985). 

There is a theoretical prediction of the power in the algebraic correction. The 
quasi-normal theories of turbulence (Kraichnan 1959; Tatsumi 1980; Lesieur 1987), 
which will be discussed in some detail in the next section, predict a = 3 in the 
far-dissipation range. This value of a is quite different from those observed in the 
numerical simulations. But it should be remembered that the dissipation range is 
restricted to k/k d < 2 ~ 3 in the simulations mentioned above so that it is not 
clear whether this discrepancy results from a failure of the quasi-normal theories 
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FIGURE 3. The energy-spectrum balance: , transfer T(fc); , dissipation 

2 uk 2 E(k)\ , departure from equilibrium \T(k ) — 2vk 2 E(k)\. The equilibrium 

between transfer and dissipation is apparent over most of the wavenumber range. 

or from the low wavenumbers considered. As a matter of fact, there is a numerical 
suggestion that the spectrum may be consistent with positive values of a less than 
3 over the wavenumber range 4 < k/k^ < 10 for a low Reynolds number ( R\ % 15) 
flow (Domaradzki 1992 private communication). 

In order to see the dominant terms of (2.13), we plot T(k) and 2uk 2 E(k) in 
figure 3. Recall that the forcing term F(k) in (2.13) is restricted to low wavenum- 
bers (k/kj < 0.1). We see that both 2vk 2 E{k) and T(k) vary exponentially with 
wavenumber and that they are in equilibrium. Their difference is less than their 
magnitude by more than two orders of magnitude over most of the wavenumber 
range (k/k^ > 1, say). 

3.2 Triad energy transfer 

The triad energy transfer function T(k\p, q) is most efficiently calculated by an 
alias-free spectral method applied to filtered fields (Domaradzki & Rogallo 1990). 
In figure 4, we plot T(k\p, q) for k/k<\ = 2 and k/k a = 3. The finest band-width 
is employed, i t e, AA: = 1. The solid and broken curves denote the positive and 
negative values, respectively. We recognize the following characteristics in T(k\p, q). 
First, there are strong dipoles at the corners q <C p & k and p q ^ k of the 
boundary. The signs of the dipoles are positive (negative) on the smaller (larger) side 
of wavenumber pj k or qfk. Second, the most significant part of T(A:|p, q) is localized 
near the boundary p + q — fc, and the thickness of this part decreases as k/k^ 
increases. The magnitude of T(fc|p, q) decreases exponentially with wavenumber 
away from this boundary. The value of T(A:|p, q) in the blank region is too low to 
draw clear curves. 

The first characteristic was also observed in the inertial range (Kida & Ohkitani 
1992) and represents the energy transferred to larger wavenumbers by nonlocal inter- 
actions. The second characteristic, on the other hand, is peculiar to the dissipation 
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p/k _ p/k 

FIGURE 4. The triad energy transfer function T(k\p,q) for the numerical tur- 
bulence. (a) k/k d = 2, (b) k/k d = 3 The solid and broken curves represent the 
positive and negative parts, respectively. The contour levels are logarithmic, rather 
than linear, and are separated by a factor of 4. 



(p - k)/k d _ ( P ~ k)/k d 

FIGURE 5. The triad energy transfer function T(k\p,q) for the numerical tur- 
bulence. (a) fc/Jfcd = 2, (b) k/k d = 3. The solid and broken curves represent the 
positive and negative parts, respectively. The contour levels are logarithmic, rather 
than linear, and are separated by a factor of 4. 


range. This arises from the rapid (exponential) decay of the energy spectrum with 
wavenumber in the dissipation range in contrast with the slow (algebraic) decay in 
the inertial range. 

Similarity in the contours evidently is not obtained over the whole domain of 
T(k\p,q) plotted in figure 4. Since, however, the dipole parts are very similar 
in figures 4(a) and (b), we enlarge the corner region and replot the contours in 
figures 5(a) and (b), respectively, with wavenumbers normalized by the Kolmogorov 
wavenumber k d instead of k. This scaling of the wavenumber is suggested by the 
closure theory (see (4.12)). 

The close resemblance of figures 5(a) and (b) implies that the shape of the dipoles 
of T( Jfc |p, q) is similar near the corners if the wavenumber is scaled by the Kolmogorov 
wavenumber. 
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FIGURE 6. The ( p/k,q/k ) domain of integration for T(k\p,q). The measure of 
scale disparity is s = max(fc,p,< 7 )/min(fc,p, q). (a) s < 2, (b) s > 2. 


3,3 Nonlocality of triad interaction 

The triad energy transfer function T(Ar|p, q) represents the energy exchange among 
three wavenumbers with magnitudes fc, p, and q . In order to express the scale 
locality of the triad interaction, we introduce the scale disparity parameter (Zhou 
1992), the ratio of the maximum to the minimum of the triad wavenumbers, 


max(fc,p, q) 
min (k,p,q) * 


(3.4) 


It follows by definition that s > 1. This scale disparity parameter measures the scale 
locality of the triad interaction. If s is smaller, the interaction is more local, and 
vice versa. In figure 6, we indicate the (fc,p, q) domains for relatively local ( $ < 2) 
and relatively non-local (s > 2) triad interactions. 

Let us denote by T(k\s)Ak the contribution to the energy transfer from triad 
interactions with scale disparity parameter between s — ^ As and s + |As. Then, 
we have 


T(k\s) = (Aib) 2 


T(k\p,q). 


(3.5) 


p»? 


In figure 7, we plot T(fc|s), obtained by summing up the terms in the r.h.s. of 
(3.5) numerically for k/k ^ =2 and kjk& = 3. It is seen that T(A:|s) may have a 
scale similarity with sk^/k. 


4. Closure Theory 

In the quasi-normal closure theories (see Tatsumi 1980; Lesieur 1987), the energy 
transfer term in the energy spectrum equation (2.13) is expressed in terms of the 
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FIGURE 7. The scale disparity of energy transfer T(k\s): , k/k& = 2; , 

k/kd = 3. (a) local scaling: s, (b) non-local scaling sk^/k. 


energy spectrum function under the assumption that the fourth-order cumulants of 
the velocity are negligible: 


7YH £ 2 ((f> . n ^(p)E(q) 

T(k\p, q ) — Bkpq { (Bkpi + Bkqp) p2 ^2 

„ E(k)E(q) „ E(k)E(p h 

-v*pq— j^2 Pp2 )' 


(4.1) 


where 

Bkpq = ( k 2 - q 2 )(p 2 -q 2 ) + k 2 p 2 (4.2) 

and 

Q 2 = 2 k 2 p 2 + 2 k 2 q 2 + 2 p 2 q 2 - k* — p* — q 4 . (4.3) 

Here Bkpq , the relaxation time of the triple moments of velocity, takes different forms 
in the various theories. In the far-dissipation range ( k,p,q k<\), however, it has 

the common expression 

Bkpq = /I 2 I 2 i in' 
u(k 2 + p 1 + q l ) 

In the far-dissipation range of statistically stationary turbulence, the first two 
terms balance in (2.13), 

T{k) = 2vk 2 E{k). (4.5) 

In this section, we omit the hat (*) because we are considering the continuous limit 
(infinite size of the periodic cube). The summation in (2.15) of the energy transfer 
function converts into the integral 


T(k) = 


I L 


T(k\p, q)dpdq, 


(4.6) 


where the integration is carried out under the condition that the three wavenumbers 
Jk, p, and q constitute a triangle. 
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FIGURE 8. The triad energy transfer function T(Ar|p, q ) for the quasi-normal closure 
theory with E(k) a (k/k < i)~ ls exp[-4.9k/k d ]. (a) k/k d = 2 and (b) k/k<i = 3. The 
solid and broken curves represent the positive and negative parts, respectively. The 
contour levels are logarithmic, rather than linear, and are separated by a factor of 
4. 


By substituting (4.6) with (4.1)-(4.4) and (3.3) into (4.5) and taking the leading 
terms in the limit of large wavenumber k > k<i, we obtain a = 3 (Tatsumi 1980). 
This value of a is not consistent with the results of the numerical simulation as 
mentioned in section 3. The reason for this discrepancy is not known at present. 
As will be discussed in the next section, other values are consistent with the theory 
if the relaxation time 9k pg is suitably modified (see (5.1)). We therefore proceed 
to examine the behavior of the triad energy transfer function expressed as (4.1) for 
the spectrum (3.3) with a = —1.6. 

In figure 8, we show the contours of T(k\p, q) expressed by (4.1) with the spectrum 
(3.3) with a = —1.6 and /? = 4.9 for both k/k a = 2 and k/k d = 3. Contrary to 
the simulation data (figure 4), we can see contour lines at very low levels clearly. 
The same characteristics of T(k\p, q) observed in figure 4 are also observed here. 
That is, (i) there are strong dipoles at the corners of the boundary of the triangle 
condition, (ii) T(k\p,q) is positive where either p or q is less than k and negative 
otherwise, (iii) the magnitude of T(fc|p, q) decreases rapidly as point (p, q) moves 
away from the boundary p + q = k y and (iv) it decreases more rapidly a s k/k 
increases. Moreover, the shape of the contours in figures 4 and 8 is very similar. 
The agreement is better for positive contours than for negative ones. 

The difference in the shape of the contours can be seen more clearly in figure 
9, which is an enlargement of figure 8 near 5 < p « It. As will be discussed in 
the next subsection, the slight difference in the shape of contours seems to be the 
main reason for the discrepancy in the behavior of T(k\s) between the numerical 
simulation and the closure theory. 

It should be mentioned here that the influence of the forcing term may not be 
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(p - k)/k d 

FIGURE 9. The triad energy transfer function T(fc|p, q ) for the quasi-normal closure 
theory with E(k) oc (A:/A:d)*" 1 ' 6 exp[-4.9fc/fcd]- (a) = 2 and (b) k/k d = 3. The 

solid and broken curves represent the positive and negative parts, respectively. The 
contour levels are logarithmic, rather than linear, and are separated by a factor of 
4. 

negligibly small. Since the fluid is forced at wavenumbers less than 3, the contours 
at q/k d < 0.1 are directly affected by the forcing term. 

So far, we have examined only the case of a = —1.6. In order to see the depen- 
dence of T(ifc|p, q ) on the value of a, we plot in figure 10 the contours of T(k\p, q ) for 
various values of a ranging from 3 to -2. It is seen that the domains of the positive 
and negative parts are insensitive to the value of a, but the shape of the contours 
other than the zero level changes depending on the value of a. For a large value of 
a, the positive and negative peaks of T(fc|p, q) are far from the corners. They move 
toward the corner as a decreases. For a > 0, the peaks are away from the corner, 
but for a < 0, a positive and a negative peak merge into a dipole at the corner (see 
(4.12)). As will be shown in the next subsection, the dominant interactions in the 
energy transfer are local for a > 1 and nonlocal for a < 1. 

4.1 Scale Disparity of Energy Transfer 

Let T(fc|s)d,s be the contribution to the energy transfer to Fourier modes at 
wavenumber k due to triad interactions for which the scale disparity parameter lies 
between s and s + d s. The contribution from all the triad interactions of scale 
disparity less than s is then written as 


J‘tWW = I J T(k\p, q)dpdq. 

1 J *k 


(4.7) 


The derivative of (4.7) gives the scale disparity of energy transfer 

T(k\s) ~ J J T(k\p,q)dpdq. 

minik.p.q) 

The integration of T(fc|s) over all s gives the energy transfer T(fc), i.e. 


(4.8) 
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p/k p/k 

FIGURE 10. The triad energy transfer function T(k\p, q) at k/kd = 1 given by the 
quasi-normal closure theory for E(k) oc (k / kd) a exp[—4.9k / k^]^ —2 < a < 3 
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The integral (4.7) is written explicitly for s < 2 as (see figure 6) 

jfW) d*'" 2 {£/'’ I!, dq+ L dp //. d9 } T(fc|p,9) (4i ° a) 


and for s > 2 as 


/•* f rp /•* r p 

I T(k\s')ds' = 2/ / dp I dq+ I dp I dq 

Jl [Jk/2 Jk-p J(,-i)k/, Jk/» 

r sk/(»-l ) rp [»k rp 1 

+ / dp dg + / dp/ dg >T(fc|p,?)(4.106) 

J k Jp/* J»k/{»- 1) J 


Differentiating (4.10a) and (4.10b) with respect to s we obtain 


(s<2) 


T(fc|s) = 2^ ^ T(k\p,k/s)dp + k T(k\sk,q)dq 

+ -^J k p^( fc b-p/ 5 )dp| 

= 2 [X f T(k\p, k/s)dp+k f T(k\sk,q)dq 

J(3-l)k/» J{»-l)k 

*u 


f9k/(9-\) 


pT(k\p,p/s)dp 


) 


(s>2) 


(4.11a) 


(4.116) 


By substituting the expression (4.1) for T(k\p,q) from the quasi-normal closure 
theory with the energy spectrum (3.3) into (4.11a) and (4.11b), we can calculate 
T(Jfc|s) explicitly. The scale disparity transfer function behaves differently depending 
on the values of a and k/k^. In figure 11(a), we plot T(fc|s) for k/k^ = 1, 2, and 
3 with a = 3. The interaction is localized around s = 2. The peak of the scale 
disparity parameter moves little as k/k^ increases. In figure 11(b), we plot T(k\s) 
for k/kd = 1, 2, and 3 with a = -1.5. The interaction becomes more nonlocal as 
k/kd increases, and the peak of the scale disparity parameter moves linearly with 
k/kd for large k/kd- 

In order to examine the wavenumber dependence of the transfer function, we 
replot it in figure 11(c) against skd/k for k/kd = 2, 3, and oo (see below for k/kd — 
oo). The scale disparity of the energy transfer seems to approach a universal form 
for large values of k/kd • The approach is faster for large values of s. 

As demonstrated in figure 10, the triad energy transfer function T(k\p,q) for 
small a (see below (4.12) for the critical value) has a double peak at the corners 
q -C p « k and p <€. q ~ k. This peak becomes steeper for larger values of k/kd 
because of the exponential decay of the energy spectrum. This enables us to make 
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3 sk^/k 

FIGURE 11. The scale disparity of energy transfer T(k\s) in the quasi-normal 
closure theory with E(k) a (k / k^) a exp[—4.9k / k^]. (a) a = 3, (b,c) a = —1.5. 
, k/k<i = 1; , k/k<i = 2; , k/k^ = 3; , k/k^ = oo. 


a local analysis around the corner to estimate the asymptotic behavior of T(k\s) 
for large values of k/k a. The triad energy transfer function (4.1) with the energy 
spectrum (3.3) behaves around q <C p ~ k as 


r (U p,,) = £ (£)”e- w ‘- (^) Yv -0>- _ j) . 

(4.12) 

The energy transfer function T(k) is calculated by integrating T(k\p, q ) with respect 
to p and q over the whole range. When a < 1, the integral is localized at the corners 
q *C p & k and p <C q % k so that the asymptotic expression (4.12) can be used. 
The result is 


/>oo /*+¥ 

T(k) = 2 / T(fc|p,g)dp 

*/0 Jk — q 

_ 2A 2 fc d / 2~ 2 

vfl a+l y(a — 2)(a — 1) 




(4.13) 
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for -2 < a < 1, where T is the Gamma function. Note that the integral is localized 
at the corners and converges only when — 2 < a < 1. 

Substituting (4.12) into (4.11b), we obtain 


T(k\s) 


A 2 k& ( k \ 

2 i /£“+ 2 VW 


e -0k/k d 


x o * a+2 * ^2<7 2 (1 — o) — |e 


+ 2<r 2 {\+o )e 2 /") ,(4.14) 


where 



(4.15) 


Note that we consider here the case k/k a ^ 1 so that 3 1 for a finite value of a. 

The asymptotic form (4.14) of T(k\s) for large values of fc/fcj is drawn for a = -1.5 
and j3 = 4.9 in figure 11(c). 


5. Discussion 

The triad interaction in the dissipation range has been investigated by analysis of 
numerical turbulence data and the quasi-normal closure theory of turbulence. The 
results of the numerical simulation suggest that the motion at the Kolmogorov dis- 
sipation scale couples directly with the smaller scales and that the triad interaction 
is nonlocal in scale, at least in the wavenumber range 0.1 < kjk& < 4 investigated 
here. The closure theory, on the other hand, suggests that the locality of the triad 
interaction depends crucially on the power a of an algebraic prefactor of the expo- 
nential decay of the energy spectrum at large wavenumbers. It is local or nonlocal 
for a > 1 or a < 1, respectively. 

In the EDQNM and related quasi-normal closure theories, the triad energy trans- 
fer function is expressed by (4.1). The form of the relaxation time Okpq differs from 
theory to theory but has the common asymptotic form (4.4) in the far-dissipation 
range where a balance between energy transfer and dissipation requires a = 3 in- 
stead of the value a & — 1.6 found in the numerical simulation. 

It is interesting, however, to note that if the relaxation time is assumed to be 
independent of the wavenumber in the far dissipation range, say equal to the Kol- 
mogorov time scale 

Qkpq <X = 71 ' 

vk& 

then any value of a (including a = —1.6) is compatible with the energy balance 
equation (4.4). As shown below, however, this is not the case in the EDQNM theory. 

In the EDQNM theory, the relaxation time takes the form 

6kpi = Kit 2 +p 2 +g 2 ) + /i(fc) + Mp) + M<z)’ (5 ‘ 2) 

where fi(k) = \(f* r 2 E(r)dr ) 1 1 2 is the eddy damping rate and A is an adjustable 
constant which may be related to the Kolmogorov constant Ck as A = 0.154Ck 3/2 
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(Andre & Lesieur 1977). Note that A = 0.37 for C*k = 1.8. The first term in the 
denominator in (5.2) represents viscous damping and the second the relaxation by 
straining motions of comparable and larger scales. The relaxation time has the 
following asymptotic forms for small and large wavenumbers for Ck = 1 -8: 


Okpq — 


2.3 

A 2 / 3 + p 2 / 3 + g 2 / 3 

1 

v(k 2 + p 2 + q 2 ) 


( k,p,q < fed) 

(*,P,9 » *d) 


(5.3a) 

(5.36) 


The integrals in (5.2) tend to the energy dissipation for large wavenumber, 


r°° 

t/2v = / r 2 E(r)dr 
Jo 


The peak of the integrand lies around r = 0.15&d> ;m d the majority of the integrand 
is covered in the wavenumber domain r < 0.5fcd (see Kida & Murakami 1987). The 
two effects are comparable at wavenumber 

k,p,q ss A 1/2 /2 ,/4 ibd « 0.5Jfc d - (5.5) 

Around these wavenumbers the relaxation time is written as 


kpq u(k 2 + p 2 + q 2 ) + 0.78 (e/u) 1 / 2 v(k 2 + p 2 + q 2 + 0.78Jt d 2 )* 


We may conclude from (5.3) and (5.6) that there is no region of constant 9k pq in 
the EDQNM theory. 

It is possible that the value of a w —1.6 observed in wavenumber range 0.1 < 
kfk<\ <4 is simply a tangent and that it approaches 3 in the limit of large wavenum- 
bers. If so, the transfer may be dominated there by local triad interactions. 
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